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The connection between f[R) theories of gravity and scalar-tensor models with a "physical" 
metric coupled to the scalar field is well known. In this work, we pursue the equivalence between 
a suitable scalar theory and a model that generalises the f{R) scenario, encompassing both a non- 
minimal scalar curvature term and a non-minimum coupling of the scalar curvature and matter. This 
equivalence allows for the calculation of the PPN parameters 13 and 7 and, eventually, a solution to 
the debate concerning the weak-field limit of f{R) theories. 
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I. INTRODUCTION 

Contemporary cosmology is faced with the outstand- 
ing challenge of understanding the existence and nature 
of the so-called dark components of the Universe: dark 
energy and dark matter. The former is required to ex- 
plain the accelerated expansion of the Universe, and ac- 
counts for about 74% of the energy content of the Uni- 
verse; the latter is hinted, for instance, by the flatten- 
ing of galactic rotation curves and cluster dynamics [J, 
and constitutes about 22% of the Universe's energy bud- 
get. Several theories have been put forward to address 
these issues, usually resorting to the introduction of new 
fields; for dark energy, the so-called "quintessence" mod- 
els consider the slow-roll down of a scalar field, thus in- 
ducing the observed accelerated expansion 0, Q. For 
dark matter, several weak-interacting particles (WIMPs) 
have been suggested, many arising from extensions to the 
Standard Model [e.g. axions, neutralinos). A scalar field 
can also account for an unified model of dark energy and 
dark matter Alternatively, one can implement this 
unification through an exotic equation of state, such as 
the generalized Chaplygin gas [5|. 

Other approaches consider that these observational 
challenges do not demand the inclusion of extra energy 
content in the Universe but, instead, they hint at an in- 
completeness of the fundamental laws and tenets of Gen- 
eral Relativity (GR). Following this line of reasoning, one 
may resort e.g. to extensions of the Friedmann equation 
to include higher order terms in the energy density p 
have been suggested (see e.g. Q and references therein). 
Another approach considers changes on the fundamen- 
tal action functional: a rather straight forward approach 
lies in replacing the linear scalar curvature term in the 
Einstein-Hilbert action by a function of the scalar curva- 
ture, f{R); alternatively, one could resort to other scalar 
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invariants of the theory Q (see [1] and references therein 
for a discussion). 

As with several other theories [§, solar system 
tests could shed some light onto the possible form and 
behaviour of these f{R) theories; amongst other consid- 
erations, this approach is based either in the more usual 
metric affine connection, or in the so-called Palatini ap- 
proach [ll| , where both the metric and the affine connec- 
tion are taken as independent variables. As an example 
of a phenomenological consequence of this extension of 
GR, it has been shown that f{R) oc i?" theories yield a 
gravitational potential which displays an increasing, re- 
pulsive contribution, besides the Newtonian term [12| . 

Another line of action lies in the comparison between 
present and future observational signatures and the pa- 
rameterized post-Newtonian (PPN) metric coefficients 
arising from this extension of GR, taken in the weak field 
limit and when the added degree of freedom may be char- 
acterized by a light scalar field [T^ . Regarding this, some 
disagreement exists in the community, some arguing that 
no significant changes are predicted at a post-Newtonian 
level (see e.g. [l^ and references therein); others defend- 
ing that f{R) theories yield the PPN parameter 7 = 1/2, 
which is clearly disallowed by the current experimental 
constraint 7 - 1 = (2.1 ± 2.3) x 10"^ 14]. This re- 
sult first arose from the equivalence of the theory with a 
scalar field model [lB|, which led to criticism from sev- 
eral fronts [3; however, a later study implied that the 
result 7=1/2 could be obtained directly from the orig- 
inal f{R) theory [l3| (see Ref. [ll| for a follow-up and 
criticism). 

Despite the significant literature on these f{R) mod- 
els, another interesting possibility has been neglected un- 
til recent times: including not only a non-minimal scalar 
curvature term in the Einstein-Hilbert Lagrangian den- 
sity, but also a non-minimum coupling between the scalar 
curvature and the matter Lagrangian density; indeed, 
these are only implicitly related in the action functional, 
since one expects that covariantly invariant terms in C 
should be constructed by contraction with the metric 
{e.g. the kinetic term of a real scalar field, g^'^ X,fiX,u)- 
In regions where the curvature is high (which, in GR, 
are related to regions of high energy density or pressure) , 
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the implications of such theory could deviate consider- 
ably from those predicted by Einstein's theory [l^. Re- 
lated proposals have been put forward previously to ad- 
dress the problem of the accelerated expansion of the 
Universe [20| and the existence of a cosmological con- 
stant Other studies have studied the behaviour of 
matter, namely changes to geodetic behaviour [19|], the 
possibility of modelling dark matter |23] , the violation of 
the highly constrained equivalence principle [1^ [3l| and 
the effect on the hydrostatic equilibrium of spherical bod- 
ies such as the Sun; also, a viability criterion for these 
generalized f{R) theories has been obtained [2^ . 

In this study, we focus on the equivalence of a theory 
displaying a non-minimal coupling of the scalar curvature 
with matter with a scalar-tensor theory; through a con- 
formal transformation of the metric, this yields a purely 
scalar theory, that is, one in which the curvature term 
appears isolated from any scalar field contribution. For 
definitiveness, we recast the theory in a form that is as 
consistent as possible with the work of Ref . [25| , and in 
close analogy with the available equivalence with f{R) 
models [isj . 

This work is divided into the following sections: first, 
we introduce the gravity model and discuss some of its 
features; then, we recast it as a scalar-tensor theory with 
a suitable dynamical identification of the scalar fields, 
and then as a scalar theory with a conformally related 
metric and redefined scalar fields. The later prompts for 
a computation of the PPN parameters /3 and 7, which is 
followed by a discussion of our results. 



II. THE MODEL 

Following the discussion of the previous section, one 
considers the following action [l9l |. 



5: 



[KhiR)+f2iR)C]^d^X 



(1) 



where k = c'^/IGttG, fi{R) (with i — 1,2) are arbitrary 
functions of the scalar curvature R, £ is the Lagrangian 
density of matter and g is the metric determinant; the 
metric signature is (— , The standard Einstein- 
Hilbert action is recovered by taking /2 = 1 and /i = 
R — 2A, and A is the cosmological constant (from now 
on, one works in a unit system where c = k = I). 

Variation with respect to the metric g^j/ yields the 
modified Einstein equations of motion, here arranged as 



-9 Sg^"" 



(3) 



The Bianchi identities, together with the identity 
(□Vy — Vu\^Fi = Rf^yV^^Fi, imply the non-(covariant) 
conservation law 



V^T^. = ^(gp./:-r^.)v^i? , (4) 

and, as expected, in the limit f2{R) — 1, one recovers the 
conservation law V^T^^, — 0. 

Since the energy-momentum tensor is not covariantly 
conserved, one concludes that the motion of matter dis- 
tribution characterized by a Lagrangian density C is non- 
geodesical. This, of course, may yield a violation of 
the Equivalence principle, if the right-hand side of Eq. 
(HI) differs between distinguishable matter distributions, 
which could be used to experimentally test the model and 
obtain constraints on A. 



A. Equivalence with scalar-tensor theory 

1 . Survey 

As in usual fi{R) models, one may rewrite the con- 
sidered mixed curvature model as a scalar-tensor theory. 
One wishes to establish an equivalent action in the Jor- 
dan frame (where the scalar curvature appears coupled 
linearly to a function of scalar fields); this may be ob- 
tained in a similar fashion to usual f{R) models (with 
f2{R) = 1), so that the equations of motion derived from 
the action functional coincide with those derived directly 
from the action Eq. ([1]). However, one cannot simply 
consider an action of the form 



(5) 



The equivalence with the action ([T]) must stem from the 
equations of motion obtained through the variation of 
with respect to the non-dynamical scalar field $; through 
a suitable definition of the functions g{^), h{^) and the 
potential V{^)^ the variation of the action ([5]) should 
yield the dynamic identification <f> = <&(-/?), so that both 
actions are equal, that is 



(FA-Fn(E B\-R.rfT ^ m g{HR))R-vmR)) = MR) , (6) 

(Fi + F2C) (^i?,. - -g,yRj - 8nGhT,y + (2) ^^^^^^^ ^ ^^^^^ 

^ [/i - {Fi + F2C) R] gf,^ + (D^^ - g^^D) (Fi + F2C) , or, conversely, 
where one defines 0^,^ = V^V,y for convenience, as well 

as Fi(i?) =d/,(i?)/di?, omitting the argument. The mat- g{^)R{^) - V(<^>) = fi{R{^)) , (7) 

ter energy-momentum tensor is, as usually, defined by ^(^) = f2{R{^)) ■ 



3 



From action ([5]), the equation of motion for $ reads 



2. Two-field scalar-tensor model 



(8) 



where the prime denotes differentiation with respect to 
$. Integrating, one gets 

5(<i>)i? + = , (9) 

which then yields an action without matter content, 



S 



(10) 



indicating that one can resort to just one scalar field 
only when deahng with pure f{R) models (that is, with 
f2{R) = 1}" This pathology has already been pointed out 
in Ref. [2J] , and in a following study |26| it is shown that 
one must resort to a second auxiliary real scalar field ip. 
In this approach, the authors chose an action of the form 



S 



[fi{4>) + f2{(b)C + i^{R-4>)]V^d^x , (11) 



so that the second scalar field acts as a Lagrange multi- 
plier, enforcing the identification (j) — R; however, varia- 
tion of the above with respect fo yields the relation 



(12) 



so that the two scalar fields are independent, if £ 7^ 
or F2 ^ (with F2 = falling back to the trivial case 
/2 = 1 and ip — Fi ((/))). The on-shell action obtained by 
replacing Eq. (fT2|) into action (fTT|) is given by 



In this study, one establishes the equivalence between 
action ([T]) and a scalar-tensor theory through Eq. (jlip . 
here rewritten as a Jordan-Brans-Dicke theory with a 
potential. 



[ijR^V{^.,ij) + f2i(f>)C]^d^x , (15) 



where (j) and ip are scalar fields and one defines 



(16) 



As discussed in the previous section, variation of the 
metric with respect to both scalar fields yields the dy- 
namical equivalence cf) = R and ip = -F'i((/') -\- F2{(j))C 
Substituting into Eq. p5|) . one recovers the action for 
the mixed curvature model, Eq. ([T|). 

With the above taken into consideration, one should 
notice that the second scalar field ip is not a function 
of the curvature (or 0) alone, but also of the matter 
Lagrangian £ (which is itself a scalar). The degrees of 
freedom of the matter fields x which appear in the La- 
grangian C{g^'\x) (including any kinetic terms) are dis- 
played in the Einstein field equations, through the corre- 
sponding energy-momentum tensor T^^: 



S 



X X 



(13) 



[/i(0) + Fi(0)(i? - 0) + [/2(<^) + F2mR - 



This differs from the unobtainable form of Eq. ([5]) in the 
dependence of F{4>,C) = J^i(0) -I- F2{(j))C on both the 
scalar field (j) and the Lagrangian matter density. 

As in usual f{R) models, a suitable conformal transfor- 
mation to the metric g^^, — F{(l), C)g^i, could be used to 
transform the above action of Eq. (fT3|) into a functional 
where the scalar curvature appears decoupled from other 
fields. However, the presence of the Lagrangian density 
in this transformation would render any comparison with 
standard scalar theories too complicated (including the 
extraction of the PPN parameters /3 and 7) ; for this rea- 
son, the second scalar field ip will be kept for the remain- 
der of this study. 

Notice also that one also could opt for the equivalence 
with a theory with just one scalar field, by taking the 
action 



RlLf 



/2(0) 



1 ^(0,V.) ^ 1 , 

1:9 + - - gf^uU) tp 



(17) 



i? and V = Fi{(l>) + F2{(j))C, 



which, introducing (I 
recovers Eq. ©. 

Using the Bianchi identities and the previous relation, 
one obtains 



1 



[v,y (0, i>) -Rv.ij- F2 (0)r^. 



Since 



(18) 



[h{^)-^Fi{^){R-<P) + f2{R)C]^d^x , 

(14) 

However, the latter is not written in the Jordan frame, 
as the scalar curvature appears in a non-linear fashion. 



we get 



V, 1/((/., V) = - Fi (0)] + (/.V, V (19) 
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I 

which, upon the substitution, ip — Fi{(f>) + F2{(f>)C and By resorting to the divergence theorem, the first integral 
4> = R, coUapses back to Eq. (|4]). may be dropped, yielding 



III. EQUIVALENCE WITH A SCALAR THEORY 

One may now perform a conformal transformation (see 
e.g. Ref. [13 )j so that the curvature appears decou- 
pled from the scalar fields 0, -0 (yielding the action in 
the so-called Einstein frame), by writing g*^ = ipg^i, = 
A~'^{ip)9fii^, with A^Tp) = tp~^^'^, one obtains 



(21) 



R = 



where □* denotes the D'Alembertian operator, defined 
from the metric 5*^. From the definition of the energy- 
momentum tensor, this implies that T*^ = A^(%li)T^^. 
Introducing the above into the action ([15]) yields 



S 



(22) 



R* - eVia* ) - 4c/ + f2A^c{A^g;,, x) 



where one defines 



U{^,i,)=A\^)V{cl,,^)lA = ^ 



<^ h{4>) 



Alp 4?^^ 



(23) 



Notice that there are two couplings between the scalar 
fields (j), Tp and matter: the explicit coupling given by 
the factor f2{(l>)A'^{ip), and an additional coupling due to 
the rewriting of the metric (in the Jordan frame) in 
terms of the new metric ^ . 

One now attempts to recast the action in terms of two 
other scalar fields, endowed with a canonical kinetic term. 
For this, one first integrates the covariant derivative term 
by parts and uses the metric compatibility relations, ob- 
taining 



f- 



6 J (^_L^yr^d4^ 



(24) 



d^x 



(25) 



One obtains the action 



S 



y/~g* d X X 



(26) 



R* 



lg*^-^^-W + hA^CiA'g;^,x) 



One may redefine the two scalar fields, so that their 
kinetic terms may be recast in the canonical way. Specif- 
ically, one aims at writing (see Ref. [25l|): 



3 V',/xV',i^ _ 3 
~ 2 



'^■^''^ - " (log ^)^^ (log V'),, ^ 2a.,,^W,. , (27) 



2 

with i, j ~ 1, 2; CTy is the metric of the two-dimensional 
space of scalar fields (field-space metric, for short), and 
<~p^ the two new scalar fields. Clearly, this kinetic 
term prompts for the identification 



1 V3, , 
V = — log ^Z- 



(28) 



and the field metric 



1 




(29) 



which indicates that only ip^ has a kinetic term. Despite 
this, it is clear that (/j^ = is a distinct degree of freedom, 
since one cannot rewrite the potential V{(p, ijj) in terms 
of one scalar field alone: 



f'^ - /i(v^)exp 



2\/3 1 



ip- X 



(30) 



2\/3 1' 



In the trivial case f2{R) = 1 or £ = 0, one gets ■0 = 
Fi(0), that is, ip^ cx \ogFi{ip'^), so that one degree of 
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freedom is lost and this potential may be written as a 
function of just one the fields. 

Since the inverse field metric fi'^ will be required to 
raise latin indexes throughout the text, one still has 
do deal with the particular form of (jy, which is non- 
invertible. In order to cope with this caveat, one is free 
to add an antisimmetric part, rewriting it as 



with inverse 



1 a 
-a 



-1 

1 a-i 



(31) 



(32) 



Clearly, no physical results can depend on the value a 
that shows in the off-diagonal part of aij - in particular, 
the values of the PPN parameters /3 and 7, as shall be 
shown. 

With this choice, the action now reads 



Alternatively, one may use the field-space metric aij 
to raise and lower latin indexes, so that a* — a^^ aj and 
= a^^B-j, rewriting Eq. jSS]) as 



and Eq. becomes 



{a^f2T* + 6^2F2C*) 6ip' , (38) 



□V ^B' - 4:TTGa'f2T* 



ATTGa'^F2C* 



(39) 



This form enables a prompt comparison with the /2 = 
1 limit, yielding the usual expression [13] for just one 
scalar field, = B - AnGaT* (with ip — ip^ and a — 
ai). 

Using the Bianchi identities, the expression for the 
non-covariant conservation of the energy-momentum ten- 
sor in the Einstein frame is also attained: 



o J2 



{g;,c* - t;,) ■ (40) 



S ^ J d'^x X (33) 

[R* - 2g*''V,,^:^^f, - AU{v\^') + f2{p')C*] . 

where one defines £* — A'^{cp^)C, with A{ip^) = 
exp (— (v^/S);/?^) . By varying the action with respect 
to the metric g*^'^, one obtains 



If, for consistency, one rewrites this in the Jordan 
frame, the conformal transformation properties of the 
contravariant derivative eliminate the first term in the 
r.h.s., and the substitution Lp^ = <j) = R yields Eq. 
Again, taking /2(i?) = 1 one recovers \'*^T*^ = aT*ip^i, 
(or, in the Jordan frame, the covariant conservation law 
VT^. = 0). 



R;.-^g;,R*^87TGf2{p^)T;^ + 



(34) 



Defining T* = g*'"'T*^ and = A-^{dA/dip'), so 
that 



ai = - — , a2 = 
one obtains, after some algebra 



(35) 



i(v^^^ = _^/,r*V + i^2/:*V , (36) 

V-ft 3 
so that the Euler-Lagrange equation for each field ip^ read 



anV = -B2 + 4:nGF2C* 



(37) 



a 6a a 



where one defines Bi = dU/dtp'^. 



IV. PARAMETERIZED POST-NEWTONIAN 
FORMALISM 

Assuming that the effect of the non-minimum coupling 
of curvature to matter is perturbative, one may write 
f2{R) = 1 + \52{R), with XS2 <C 1, so that the current 
bounds on the equivalence principle are respected. Sub- 
stituting into Eq. one gets, at zeroth-order in A, 



(41) 



which amounts to ignoring the f2{p'^) factor in the action 
([55)1 , so that /2 manifests itself only through the coupling 
A^{ip^) present in T* , and the derivative of ip^ (since 
p^ oc logip and -0 = Fi -|- F2C). In this case, the matter 
action reads 



A\^)C{A^{i^)g*x)V^d'^ 



(42) 



If there is no characteristic length ruling the added 
gravitational interaction (that is, both scalar fields are 
light, leading to long range interactions), this manifesta- 
tion of a "physical metric" in the matter action allows 
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one to resort to calculate the PPN parameters /? and 7 



(43) 

where aj^i = daj/dip^ and = a^a' = a''^aiaj; the sub- 
script indicates that the quantities should be evaluated 
at their asymptotic values ipQ , related with the cosmologi- 
cal values of the curvature and matter Lagrangian density 
(as shall be further discussed in the following section). 

Since the a^'s are constant for z = 1, 2, one gets that 
aj^i = 0; hence, the PPN parameter /3 is unitary. More- 
over, since ai = (ai,0) is orthogonal to a* = (0, ai/a), 
one gets = 0. Hence, the PPN parameter 7 is also uni- 
tary. As argued previously, these results do not reflect 
the particular value of a chosen in the antisymmetric part 
of aij. 

In light of the overall discussion presented, it is clear 
that the two degrees of freedom embodied in the two 
independent scalar fields Lp^ and (p'^ stem not only from 
the non-minimal curvature term in action ([T]) , but also of 
the non-minimal coupling of R to the matter Lagrangian 
density. Given this, it is clear that a "natural" choice 
for the two scalar fields (in the Jordan frame) would be 
(j) — R and ip = C: this more physical interpretation 
ab initio comes at the cost of more evolved calculations, 
since both redefined scalar fields in the Einstein frame 
will depend on (j) and ip. This less pedagogical approach 
is deferred to Appendix [K\ 



more elaborate calculations would have to be performed 
in order to extract the PPN parameters (3 and 7. 

It is important to highlight that a naive analysis of the 
model under scrutiny might predict no difference between 
the PPN parameters arising from the trivial f2{R) — 1 
and the non- minimal f2{R) 7^ 1 cases: indeed, since 
this function is coupled to the matter Lagrangian den- 
sity one might expect that, outside of the matter distri- 
bution (£ = 0), the theory would collapse back to the 
usual f{R) scenario. However, this misinterpretation is 
resolved by Eq. which is evaluated at the cosmo- 

logical values (f}^: for this reason, the relevant value for 
C is given by the Lagrangian density of the overall cos- 
mological fluid, not of the local environment. The issue 
of suitably identificating this contribution is discussed in 
Ref. m. 

Finally, note that the result /3 = 7 ~ 1 is approximate, 
since it corresponds to dropping the term /2(v3^) in action 
P3p: a future work should consider a pcrturbative ap- 
proach, which could perhaps yield a small, A-dependent 
deviation from unity, thus marking a clear (even if small) 
departure between the model studied here and General 
Relativity. 

APPENDIX A: ALTERNATIVE FORMULATION 
OF THE TWO-FIELD EQUIVALENT 
SCALAR-TENSOR THEORY 

We present here an alternative formulation where one 
chooses instead to express action ^ through the equiv- 
alent expression 



V. DISCUSSION AND CONCLUSIONS 

The result 7 = 1 is key to our study: it is clear that, in 
the standard f2{R) = 1 theories, does not vanish (it 
is a purely algebraic, not matricial result and = 1/3), 
and the resulting PPN parameter 7 = 1/2, which vio- 
lates well-known observational bounds! A more thorough 
discussion on the ongoing debate concerning the value of 
the PPN parameter 7 for /(i?) theories is deferred to Ap- 
pendix |B]- with special focus to what is believed to be 
a misconception in the identification of the equivalence 
with a scalar-tensor theory. 

In the /2(i?) ^ 1 case, the added degree of freedom 
that a non-minimal coupling of curvature to matter im- 
plies yields not one, but two scalar fields: as a result, a 
two-dimensional field-space metric aij arises; from the re- 
definition of the fields necessary to absorb non-canonical 
kinetic terms after the conformal transformation to the 
Einstein frame, it follows that this enables a vanishing 
= aiUjU^^ term, yielding no post-Newtonian obser- 
vational signature that discriminates these models from 
General Relativity. However, this conclusion is valid only 
in zeroth-order in A: if more terms are allowed, the non- 
covariant conservation law for the energy-momentum 
tensor is no longer of the form treated in Ref. i25i|, and 



S = j [F{^,^)R-V{^,^) + f2{cl>)C{g^.,x)]V^d''x , 

(Al) 

where one defines 



FicP,^) = F^icP) + F2{(l>)^ , 

F(0,V) = 0n0»~/i(0) • (A2) 

One obtains the equivalence with the model under 
scrutiny by writing the equations of motion for the scalar 
fields, 

F2{<l>) (£ - ^) + [Fi{^) + F^{c^)i:] (i? - 0) = , (A3) 

F2mR-c^) = , 

implying that cf) — R (or i^2(0) = 0^/2 = 1, the trivial 
result) and, therefore, ip = C . Substituting into Eq. 
(jAip . one recovers the action for the mixed curvature 
model, Eq. (P). 

It should be stressed that the GR limit fi = R and 
/2 = 1 disables the identification (f) — R: one may write 
/i(0) =(j) + eSi{(j)) and /2(0) = 1 + AJ2(0), so that Eqs. 
(|A3P become 
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XS'M (£ - ^) + [eS'licj,) + XS'^{^)4>] (i? - ^) = ,(A4) 

XS'{(f>){R-(j)) = , 

and taking the limit — > gives a trivial identity. For 
this reason, it is misleading to insert the above approx- 
imations in results that stem from the scalar field ap- 
proach to f{R) theories, and only then consider the GR 
limit: the formalism itself breaks down at its inception. 
For this reason, one concludes that one cannot simply 
take the limit e ^ and argue that, as a f{R) theory 
collapses back to GR, so should the PPN parameter 7 
approach unity, which does not happen if 7 = 1/2 and 
does not show a dependence on e. 

Variation of the action (jAip with respect to the metric 
yield the Einstein equations. 



which, introducing (j) = R and if) = C, recovers Eq. ([2]). 

Using the Bianchi identities and the previous relation, 
one obtains 

V^T^, = (A6) 

12(0) 

Since 

V,V{^,^)= (A7) 

and 

V,F(0,^)= (A8) 
(Fi'(0)+F^(0)V^)V,</) + ^^2(</>)V.V , 



F{c^,i^) [ R^, - -g^^R ) = 8iiGh{c^)T^u (A5) 



we get 



v^r„ 



/2W 



(A9) 



r 



which, upon the substitution tjj ~ £ and (j) — R, collapses 
back into Eq. (g]). 



1. Equivalence with a scalar theory 

In this alternative formulation, the adequate confor- 
mal transformation to decouple the scalar curvature from 
the scalar fields 4>, tp is given by g*^ — F{(j),i{})gfj,^ = 
A-'^{(j),ip)gf^^, with A{(l),ip) = V")- One obtains 



-9 = F- 
R = F 



(AlO) 



R* - 6%/Fn* 



Introducing the above into the action (jAip yields 



V— 5* d X X 



(All) 



R* - 6%/Fn* 



iU + hA^C{A'g*,x) 



where one defines J7(0,'0) — A'^ {(j), 'ip)V {(j), iIj) / 4:. 

One now attempts to recast the action in terms of 
two other scalar fields, endowed with canonical kinetic 



term. For this, one repeats the integration of the covari- 
ant derivative term by parts and uses the metric compat- 
ibility relations, finally obtaining the action 



S = 
R* 



'—g* d X X 



(A12) 



^g*,^E^^4U + hA^C{A'g;^,x) 



As one aims to allow for an immediate comparison with 
the /2 = 1 scenario, it is interesting to isolate the con- 
tributions arising from the non-minimal scalar curvature 
term and from its coupling with matter as clearly as pos- 
sible, namely 



tPF2 



log F = log (Fi + F2V') = log Fi + log ( 1 
assuming that Fi 7^ 0. If one defines 



(A13) 



^-—logF.W , ^=—log[l + ^r^ 

(A14) 

then the comparison is transparent: when /2 = 1, the 
second scalar field vanishes, and the first scalar field co- 
incides with the usual redefined scalar field (p^ = log Fi 
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arising in fi{R) models (although in many studies some 
terms are overlooked, see Appendix |B]). Also, when 
/i — R, the first scalar field ip^ vanishes. 

For the particular choice of fields Eqs. jAMj, one 
obtains cross-products between the derivatives of ip^ 
and so that the aij field-space metric displays non- 
vanishing off-diagonal elements; since one has absorved 
the numerical factors in the rcdifinition, this field-space 
metric is trivially given by cr.y = 1 for i, j = 1,2. How- 
ever, this field-space metric is not invertible, a problem 
already found in the main body of this work; as before, 
this issue may be surpassed by adding an antisymmetric 
part to this field metric, adopting instead the form 



1 1 
1 1 



l-a 
a - 1 



1 2-a 
a 1 



with inverse 



1 a- 2 
-a 1 



(A15) 



(A16) 



As already discussed, the physical results will not depend 
on the value a. 

With this choice, the action now reads 



and, given definitions Eqs. (jA14|) . 



^^^^=-2aS. (A21) 



difi^ 3 F{ 
The Euler-Lagrange equation for each field reads 



2Fi Fo 
f2T* - -^5aC 
^1 



(A22) 



In the /2 = 1 limit, one recovers \Z\*'P — B — AnGaT*, as 
before. 

Using the Bianchi identities, the expression for the 
non-covariant conservation of the energy-momentum ten- 
sor is 



^1 h 



(A23) 

As discussed regarding Eq. ([40]) . the first term of the 
r.h.s is elliminated in the inverse conformal transforma- 
tion back to the Jordan frame, and one obtains an ex- 
pression that collapses to Eq. (|4|), after a suitable ma- 
nipulation. 



•5 - y d'^x X (A17) 

where one defines C* = A'^{ip^ , ip'^)£ as before, but using 
the new redefined fields tp*; notice that the non-minimal 
coupling is written in terms of ip^ alone. The Einstein 
field equations are 



(A18) 



and variation of the matter action with respect to each 
scalar field yields, after some algebra 



2Fi F2 

f2T* - 

^1 



where, as before, one defines 



(A19) 



Sip' = 



dlogA IdlogF ^/3_ 

"^ = ^^^-2^;^^-- = " ■ ^^^^^ 



2. Parameterized Post-Newtonian formalism 

Once again, assuming a perturbative effect of ,f2{R) 
1 + A(52(i?), with A(52 ^ 1 yields, to zeroth-order in A 



^*^T*^a,T*pl, , 



(A24) 



which amounts to ignoring the f2{<P) factor in the ac- 
tion (IA17p . so that /2 manifests itself only through the 
coupling A^{(t),'ip)] in this case, one may write 



(A25) 



The computation of the PPN parameters (5 and 7 is 
very similar to the one already performed: since a^'s are 
constant for i = 1, 2, one gets that aj,i = and /3 = 1. 
One also obtains 



(A26) 



(1 



[1 1] 



1 a-2' 




' 1 ' 


-a 1 




1 



= 



independently of the value a, as already discussed. 
Hence, one confirms that j — 1. As required, the partic- 
ular choice of field-metric a and scalar fields (initially in 
the Jordan frame) and '0 do not change the obtained 
results /3 = 7 = 1. 
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APPENDIX B: DISCUSSION OF PREVIOUS 
RESULTS 

In this appendix one addresses the issue of the f2{R) — 
1, = f{R) case, very much discussed in the liter- 

ature (see e.g. Ref. [7|) - in particular, the claim that 
f{R) models characterized by a light scalar field still al- 
low for the PPN parameters (3 and 7 to be close to unity. 
The action for this model is, in the Jordan frame, given 
by 



S = J [f{R) + C{g,,.,x)]V^d'x , (Bl) 
which is dynamically equivalent to 

5 - y - + C{g,,,, x)] d^x , (B2) 

with the usual identification (j) = R, and the definition 

y(0) = 0F(0)-/(^). 

A conformal transformation g*^^ — Fg^^, — A ^g^v, 
with A{4>) ~ F^^/^ ((/)), yields the action, in the Einstein 
frame. 



S 



R* 



^—g* d X X 

^-g*'^'^^-4U + A'C{A'g%,x) 



(B3) 



where U{(f>) = A*{(l>)V{(p)/A. 

One may redefine the scalar field as 

ifi^^log - -V3 log A , (B4) 
so that the action becomes 



J d'^x X (B5) 

[R* - 2.g*''>,^^,, -AU + A^C{A^g;,, x)] , 

and the matter action depends not on the Einstein metric 
(7* J,, but on the original Jordan metric = A^g*^^: 

Sm ^ J A^CiA'g;,,x)V^d'x^ (B6) 



i^{9tiu,x)\^d'^x 



One obtains 



dlog A \d log F \/3 



dip 



2 Olp 



(B7) 



identical to the previous result for aa- However, in this 
case the field-space metric a is one dimensional, and sim- 
ply given by an = a^^ — 1, so that ~ 1/3. One 
obtains a.j^ — 0, so that 



7-1 = -2 



a a , 



1 + 



0^/3 = 1 , (B8) 
1 1 



7 = 



2 ' 



This indicates that general scalar-tensor theories with 
no a priori kinetic term for the long range scalar field (in 
the Jordan frame) are incompatible with observations. 
As the above example shows, the equivalence between 
f{R) theories and such models falls within this category, 
and is therefore observationally ruled out. Furthermore, 
the result = 1/3 enables, for the case of a single scalar 
field (see Ref. [25|), the identification of the Brans-Dicke 
coupling parameter 



2cj + 3 = 



w = 



(B9) 



which shows that f{R) models may also be recast as a 
(sometimes used) generalized Jordan-Brans-Dicke model 
with no kinetic term, 



5* 



[^R - 1/(0) + C] 



-gct^x 



(BIO) 



with the dynamical identification = F{R) and a suit- 
able potential ^(0) ^ R{(j))F{R{(j))) - /(-R((?!))). 

In several papers in the literature this equivalence with 
a scalar-tensor theory is given by the action 



S = J d^x X (BU) 

[F{cP)R ~ Z{^)gf^-'^^^^^, - y(0) + £(5^., x)] , 

with Z{(j)) = 1; for later convenience, one retains the 
kinetic function Z(0). As discussed above, one can opt 
by an equivalent Jordan-Brans-Dicke theory with a scalar 
field dynamically identified through = F(R), and no 
kinetic term, that is, u = 0. 

It is easy to verify that variation of the action with 
respect to the scalar field 4> will yield terms involving 
Z'{4>) (which vanishes, in the usual approach Z((j)) = 1 
and the four-dimensional D'Alembertian operator, sim- 
ilarly to the classical Klein-Gordon equation. For this 
reason, the presence of a kinetic term in the above action 
implies that the dynamical identification (j) = R (arising 
from the equation of motion of the scalar field (f) fails. 

Moreover, the above conformal transformation g*^^, = 
F{<l))9t^i^ yields 
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5" = 



11 

2^ p2 



-g* d X , 



(B12) 



using the previous result relating R* and R, as well as 
gMi' ^ F{(j))g*t"' and y=g = F-^{(j))y/^. The usual 
redefinition of the scalar field follows (erroneously) , 



d^Y 3 fd log F{(j))V Z{(t>) 



d(j) J 4V d(t) J 2F{(j)) 
The usual redefinition J^] is often presented as 



difiY 3 fd log F{(f>)y 1 



(B13) 



(B14) 



which clearly corresponds to Z((p) = 1; this yields the 
canonical action 



7-1 

(3-1 



Z{^)F(4>) + 2F'(0)2 
1 F{^)F'{^) 



A2Z{(j>)F{<j)) + 3F'{(t>)^ d(j3 



(B17) 



which, for Z{(j>) = 1, yields the PPN parameters /3 and 7 
used in Refs. p^ . 

Hence, it appears that the PPN coefficients calculated 
for a wide variety of f{R) models, and obtained by the 
dynamical identification cj) = R, are inaccurate: by rein- 
stating the correct factor Z{(p) = into Eq. (|B13[) . one 
recovers 



S = J d a; X (B15) 

[R* - 25*^>,^^,, - 4U{ip) + A^i^)C{A^iip)g;,, x)] , 

which can be matched with action Eq. (IB12P through 
the relation 



dip 
90 



(B16) 



d log F{(l)) 
d4> 



However, the above shows that a proper treatment 
should use Z{(f)) — 0, as there is no intrinsic kinetic term 
in the original, Jordan frame theory Eq. (jBlip . This 
redefinition of the scalar field will affect the calculation 
of a = 9 log A/ dip and, as a consequence, yield incorrect 
predictions for the PPN parameters (3 and 7; in particu- 
lar, one obtains a dependence on F{(j)) which would oth- 
erwise be missing. This can be seen from the following 
expressions [29j . 



dip 



3 (dJogFi^y 
A\ d4> ) 

V3, 
2 



(B18) 



log F[cj>) 



so that the calculations for the PPN parameters 7 = 
1/2 and (3=1 follow as argued previously - and Eqs. 
(jB17[) clearly show. Moreover, notice that (as already 
discussed) one cannot simply identity GR with the limit 
e ^ of a model with f{R) = R + edi{R), since the 
corresponding limit of Eqs. (jB17p (taking the correct 
factor Z{(j}) = 0) is ill-defined for F'(0) = eS"{(j}) 0. 
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